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We show that the leading-power B meson wave function can be extracted reliably from the photon 
energy spectrum of the B — > jlv decay up to 0(l/ml) and 0(o? s ) uncertainty, mi, being the 6 quark 
V) . mass and a s the strong coupling constant. The 0(1/ nib) corrections from heavy-quark expansion 

can be absorbed into a redefined leading-power B meson wave function. The two-parton 0(l/rrib) 
corrections cancel exactly, and the three-parton B meson wave functions turn out to contribute 
at 0(l/ml). The constructive long-distance contribution through the B — » V — > 7 transition, V 
being a vector meson, almost cancels the destructive 0(a s ) radiative correction. Using models of 
the leading-power B meson wave function available in the literature, we obtain the photon energy 
spectrum in the perturbative QCD framework, which is then compared with those from other 
approaches. 
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The two-parton leading-power (LP) B meson wave function (distribution amplitude) ci+ plays an essential role in a 
perturbative analysis of exclusive B meson decays based on fcy factorization theorem 0, 0, |]| (collinear factorization 
theorem 0, IE IE 13 ) Its behavior certainly matters, and has been investigated in various approaches recently. 
Models of the distribution amplitude 4>+(x) with an exponential tail in the large x region have been proposed [Tflj, 
O 1 where x is the longitudinal momentum fraction carried by the li ght spectator quark. Neglecting three-parton dis- 
tribution amplitudes in a study by means of equations of motion [ill Il2| , (x) was found to be proportional to a 
step function with a sharp drop at large x fH| . The wave function (f>+(x, kx), where kx is the transverse momentum 
carried by the light spectator quark, was also derived in the same framework [f3j . All these models depend on at 
least one shape parameter, whose determination requires experimental inputs from exclusive B meson decays. 
. £h ' I n t nis paper we shall show that the radiative decay B — > ylv provides the cleanest information of the LP B meson 
wave function </>+. This mode has been widely studied in[l|l[3[ll[Ii[ll[ll[3 

due to different motivations: for extracting the B meson decay constant fs and the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element \V u b\, for demonstrating the next-to-leading-order (NLO) calculation and the proof of QCD 
factorization theorem, for deriving resummation of large logarithmic corrections, for studying long-distance effect and 
annihilation mechanism,.... The subject on the extraction of the B meson wave function from the B — > jlv data 
has not yet been discussed. It will be shown that two-parton next-to-leading-power (NLP) [0(1/ mb")] corrections 
cancel exactly, rrif, being the b quark mass. The contributions from higher Fock states, the three-parton B meson 
wave functions, turn out to be of 0(l/m|). The constructive long-distance contribution through the B — » V — » 7 
transition, V being a vector meson, almost cancels the destructive 0(a s ) radiative correction, a s being the strong 
coupling constant. The effect from bremsstrahlung photon emissions vanishes like the lepton mass because of helicity 
suppression. Therefore, the extraction of </>+ from the measured photon energy spectrum of the B — > jlv decay suffers 
only 0(1 /ml) and 0(a 2 s ) uncertainty. 

We identify and discuss the higher-power corrections to the B — > jlv decay in Sec. II, and calculate the long- 
and short-distance effects in Sec. III. Section IV is the conclusion. The hard kernel associated with the three- 
parton distribution amplitudes is derived in the Appendix, whose explicit expression is necessary for demonstrating 
the smallness of the higher-Fock-state contribution. Our conclusion differs from that drawn in j29), in which the 
semileptonic decay B — > -kIv was regarded as a more ideal process for extracting the B meson wave function. The 
argument is that the radiative decay B — * jlv, receiving a large long-distance uncertainty, does not serve the purpose. 
As stated above, this long-distance effect is in fact cancelled by the 0(a s ) short-distance one almost exactly. 
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FIG. 1: Lowest-order diagrams for the B — > jlu decay. 



II. HIGHER-POWER CORRECTIONS 



In this section we identify and discuss higher-power corrections to the B — ► decay. The B meson momentum 
Pi and the photon momentum P2 are parameterized, in the light-cone coordinates, as 

P 1 = ^|(1,1,0 T ), P 2 = ^(0,77,O T ), (1) 

respectively, where r\ = 2P 7 /mB, mB being the B meson mass, denotes the photon energy fraction. The decay 
amplitude is decomposed into 

- e (7 (P 2 , ct) |«7m(1 " 7s)&| P (Pi)) = €^ap^v a Pi F v (q 2 ) + i [e* Ttl (v ■ P 2 ) - (e* T ■ v) P 2fl ] F A (q 2 ) , (2) 

where e is the electron charge, e-r the polarization vector of the photon, v = P\/m,B the B meson velocity, and 
q 2 = (Pi — P2) 2 = (1 — v) m B the lcpton-pair invariant mass. The decay spectrum is then given, in terms of the form 
factors Fy,A, by 

f = "i^*' 2 "^ 1 - [Fl(q 2 ) + F 2 (q 2 )] , (3) 

with a = e 2 /{Ait) and the Fermi constant Gf- 

The collinear factorization theorem for the form factors Fy,A m the large 77 region has been proved in [2*3. 127) . which 
are expressed as the convolution of hard kernels with the B meson distribution amplitudes in the momentum fractions 
x of the light spectator quark. A hard kernel, being infrared-finite, is calculable in perturbation theory. The B meson 
distribution amplitudes, collecting the soft dynamics in exclusive B meson decays, are not calculable but universal. 
In the framework of factorization theorem, there are four sources of higher-power corrections to the B — > jlv decay: 

1. The heavy-quark expansion of the heavy- light current in Eq. @, 

«7/i(l - 75)6 -> U7 M (1 - 7s)/i + ^-"7m(! " 75)* ph + 0{l/m 2 ) , (4) 

where the operator D represents the covariant derivative, and the rescaled b quark field h is related to the full 
field b by 

h(z) = l -±l e —b{z) . (5) 

The factorization of the transition matrix element associated with the first (second) term in the above expansion 
leads to the LP (NLP) B meson distribution amplitudes. 

2. The higher-power interactions in the Lagrangian of the heavy-quark effective theory (HQET). The insertion of 
the HQET interactions, 



0x = ^-h{iD) 2 h , 2 = , (6) 

771b 2?7Jfc 



into the transition matrix element associated with the first term in Eq. Q yields 0(l/m},) corrections. We 
mention that there exists an alternative heavy-quark effective theory, in which the higher-power corrections are 
formulated in a different way |30| . 
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3. The higher Fock states of the B meson. The nonlocal matrix element, 

{Q\u{z) 9 G a p{uz)h{Q)\B(Pi)) , (7) 

defines the three-parton distribution amplitudes, where G a p(uz) is the gluon field strength evaluated at the 
coordinate uz, < u < 1. The additional valence gluon, attaching internal off-shell quark lines, introduces one 
more hard propagator, i.e., one more power of l/nib- 

4. The subleading parton-level diagrams (hard kernels). The two-parton lowest-order hard kernels are displayed 
in Fig. 1, where the upper quark line represents a b quark. It is easy to observe that Fig.^a) [(b)] represents 
the LP (NLP) hard kernel, since the internal quark line is off-shell by m^A (m|) with A being a hadronic scale, 
such as the mass difference raj — mj. 



A. Heavy-quark Expansion 

The factorization of soft dynamics from the transition matrix element associated with the first term on the right- 
hand side of Eq. (gj, 

(j(P 2 ,e T )\u^(l- l5 )h\B(P 1 )) , (8) 
leads to the nonlocal matrix element [Tof. 

^^e^ + ^- k --)<0|^(^)^(0)|B(P 1 ))=i^{( ! P 1+ mB)7 5 [^ + $ + (A ; )+ . (9) 

which defines the two-parton LP B meson wave functions <I>±, with the null vectors n + = (1, 0, 0t) and n_ = (0, 1, 0t), 
and the light quark momentum k. Because the photon momentum P 2 has been chosen in the minus direction, the hard 
kernels for the form factors Fy,A are independent of the component k~ , which becomes irrelevant. We construct the 
B meson distribution amplitudes (j)±(x), x = k + / P^ , from the B meson wave functions 4>±(x, kr) = P^^±(xP^ , kx) 
by integrating the latter over kx, 

cf> ± (x) = J Skr^x.kT) . (10) 

The dependence of 4>±(x) and of <fi±(x, fc*r) on the renormalization scale fj, has been suppressed. 
Define the moments of the B meson distribution amplitude </>+(x), 

A ee / dx^^- , Ai = / dxcj> + (x) . (11) 



The asymptotic behavior of </>+ (x) has been extracted from a renormalization-group equation, which exhibits a decrease 
slower than 1/x That is, the normalization Ai of the B meson distribution amplitude is divergent after taking 

into account the evolution effect. It has been argued that a non-normalizable B meson distribution amplitude does 
not cause a trouble in practice :33], since only the inverse moment Ao matters at LP |2^,0|, which is convergent. 
Note that a hard kernel would not be as simple as 1/x at higher orders in a s , and information of more moments is 
also necessary. In the following discussion we shall neglect the evolution effect, and assume that 4>+(x) is normalized 
to unity, i.e., Ai — 1. Since the B meson distribution amplitudes absorb soft dynamics, the light quark momentum k 
is of 0(A). We then have the relative importance Ai/Ao ~ A/rrib for x ~ 0(A/mfe). 

The factorization of soft dynamics from the transition matrix element associated with the second term on the 
right-hand side of Eq. gives the nonlocal matrix element, 

(0\u p (z)i ph s (0)\B(P 1 )) . (12) 

The factorization of the transition matrix elements with the insertion of the 0(l/m&) interactions in Eq. © into 
Eq. (JHJ) leads to, 



(0\i J d 4 yT[^^)/i 5 (0)O 1 , 2 ( 2 /)]|5(P 1 )) • (13) 

The contributions from Eqs. (|12|) and l|13l) can be absorbed into the nonlocal matrix element, 

(0|w p (z)6 5 (0)| J B(P 1 )) , (14) 
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FIG. 2: Three-parton contribution to the B — > decay. 

where the rescaled 6 quark field h has been replaced by the full field b. It is easy to check that the heavy-quark 
expansion of Eq. (|14|l generates Eqs. Ijl2(l and jS) . This absorption makes sense, because Eqs. Ijl2|l and 
concerning only the initial 6 quark, are universal for all exclusive B meson decays. The decomposition in Eq. 
still holds, but the B meson distribution amplitudes <I>±, redefined by Eq. (|14|l in terms of the full field b, exhibit a 
renormalization-group evolution different from that in Eq. @ |35| . 

B. Three-parton Distribution Amplitudes 

We explain that the nonlocal matrix element in Eq. is negligible in the current accuracy: the three-parton 
distribution amplitudes, whose contributions to the form factors are supposed to be of 0(l/m;,), turn out to appear 
at l/m\. The three-parton distribution amplitudes $y, $a, Xa, and Ya in coordinate space are defined via the 
decomposition, 

{O\u p (z)gG a0 (uz)n'Lh 6 (p)\B(I\)) 



-itJapn^&viti u) —n_ a X A (t, u) H — ^_ Ya (t, u) 



(15) 

Sp 



with the variable t = v ■ z. The corresponding hard kernels arise from the contraction of all the structures T — v a fi-., 
v ■ n_7„,..., in Eq. (|15|) with Fig. El written as 

^(3) HfAu fi+l a {P 2- h-u h) ~ u{?2- h-u h)l a jl+h^q - 7s)( A + mghsT} 

a [(P 2 - h - uk 2 ) 2 } 2 ' 1 j 

where k\ (ki) is the momentum carried by the light quark (gluon). The derivation of the above expression is referred 
to the Appendix. 

For r = v ■ n_7 Q , Eq. vanishes because of • n + — e* T ■ (Pi — k\ — uki) = 0. Express a a /3n^_ = i(n- a — fi-^ a ), 
in which the first term has the same structure as of Xa- The second term fi-ja renders Eq. vanish for the same 
reason. For the other structures v a fi-, n- ai and n_ Q we always have 7° = 7 + . Once 7" = 7 + , Eq. (|16(l is 
proportional to 

(3 ) ^ P^jh+Uk,) 

a [(P 2 • (fci + ufc 2 )] 2 ' 1 J 

Note that fc^ and k£ are of O(A), and that the moments of the three-parton B meson distribution amplitudes are at 
most of 0(A 2 ) Therefore, when convoluting Eq. (|17|l with the three-parton distribution amplitudes, the resultant 
contribution to the form factors Fva is of 0(A 2 /m 2 ) compared to the LP one from Fig.^a). The same higher Fock 
state has been shown to giv e a power-suppressed correction to the B — > jlv decay in the framework of soft-collincar 
effective theory (SCET) [23. With a similar reasoning, the three-parton B meson wave functions also contribute at 
0(l/m 2 ) in kx factorization theorems. We emphasize that the three-parton B meson wave functions are relevant in 
the NLP analysis of the B — > ir transition form factors. This is the reason the B — > 7Z1/ decay is a cleaner mode than 
the B — > ttIv decay for determining the LP B meson wave function. 
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C. NLP Hard Kernels 



Contracting Fig. ^ with the two structures in Eq. ©, we get the quark-level amplitudes, 



m; 

Ml 



i tr[f{P 2 - jf) 7/1 (l - 75)(A + m B ) l5 jC+] 



4V2 
i ^[7^(1-75) 



{Pi - ky 



4V2 



(q — k) 2 — m b 



(18) 



in M^ b being replaced by n_. As stated above, Fig. QJa) is LP, because of 



and M. a b with the null vector 
(P 2 - kf = -2P 2 ■ k ~ 0(m b A), and Fig. [2b) is NLP, because of (q - k) 2 - m\ = -2P 1 ■ P 2 ~ 0(m 2 ). The 
contribution from Fig. ^b) has not yet been considered in the literature. We shall neglect the mass difference 
between the B meson and the b quark in M. b '~ in our analysis accurate up to NLP. 
The collinear factorization formulas for Fy a are written as 



F v{ A)(q 2 ) = h 



<j) + (x)H+ {A) (x, 77) + 0_ {x)H v(A) (x, 77) 



(19) 



where the hard kernels H are extracted according to Eq. (J2J by keeping only the longitudinal component k^ 
Eq. l(T%|) . In terms of the LP and NLP moments in Eq. (fTT]> , Eq. (|T^|l becomes 



F v .a(i 2 ) = 



h 



r\m B 



A ± ( 1 + - ] Ai 



(20) 



in which the coefficient 1 of Ai comes from Fig.^a) and I/77 from Fig.QJb). It has been mentioned that the equality 
of Fy and Fa at LP is attributed to the spin symmetry in the large- recoil region [2(| • The coefficient 1/rj implies the 
increase of the subleading-power correction with the decrease of the photon energy. This is why a perturbation theory 
is reliable only in the large 77 region. The distribution amplitude </>_(#), contributing only through the normalization 
of the combination, 



dx[4>+(x) — (j)-(x)] = , 



(21) 



disappears from Eq. H2U|) . As shown in Eq. (|20|l . the normalization Ai does appear at NLP, which is divergent under 
the evolution. This is another example that the QCD-improved factorization (QCDF) approach based on collinear 
factorization theorem breaks down at NLP |34l |3(j . 
The decay spectrum in Eq. becomes 



dT 



*G% \Vub 
48tt 2 



-f B rnl(l 



77)77 



A 



77 



(22) 



The above expression indicates that the NLP terms for the spectrum have cancelled, and only the 0(\/m 2 ) term A^ 
is left. In this case we can estimate the 0{l/m 2 ) effect using the models for the B meson distribution amplitudes 



available in the literature 13. 3 



4>±{x) 



\±(x - X) 
2A 5 



6(x)6{2\ - x) , 



(23) 



with the shape parameter A = A/to;,. The value of A has been found to range between 0.5 and 0.7 GeV |25l 1381 l39j. 
which corresponds to A = 0.1 ~ 0.15 approximately. Certainly, there are other models of the B meson distribution 
amplitudes (see 0D- 

Employing the inputs a = 1/137, G F = 1.16639 x 10~ 5 GeV" 2 , \V ub \ = 3.9 x 10~ 3 , f B = 190 MeV, and m B = 5.28 
GeV, we derive the photon energy spectra of the B — > -flu decay for A = 0.1 and for A = 0.15 in Fig. EI The specific 
models in Eq. I|23|l lead to the relation A1/A0 = A. Therefore, the subleading-power term is indeed negligible at 
large 77, whose contribution is around 5%. However, this term diverges quickly at small 77, breaking the perturbative 
expansion in 1/mb- The form factors Fv,a in Eq. (|20|l contain a dominant monopole component proportional to Ao/77, 
and a small dipole component proportional to A1/77 2 , which is important only at small 77. This is the reason one 
always obtains a symmetric spectrum in 77 at LP from a perturbation theory |20j as shown in Fig. |3J To generate an 
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FIG. 3: Spectra Spectra in units of GeV from collinear factorization with the solid (dashed-dotted)line corresponding to 
the LP contribution for A = 0.1 (A = 0.15), and the dashed (dotted)line to the inclusion of the NLP contribution for A = 0.1 
(A = 0.15). 



asymmetric spectrum, the dipole component must be enhanced as postulated in |16L l26|. Therefore, an asymmetric 
spectrum signals an important NLP contribution, ie., a breakdown of factorization theorem. 

It has been explained that the undesirable feature of the B meson distribution amplitude under evolution is a 
consequence of collinear factorization, which can be removed in hr factorization |35| . The evolution effect on the 
fcr-dependent B meson wave function was also studied in |41j. Moreover, applying kx factorization theorem to the 
B — > jlv decay, which has been proved in |3( , we can extend the spectrum to lower rj as demonstrated below. Keeping 
both the longitudinal momentum k + and the transverse momentum kr in Eq. I|18|) , the hard kernels in Ict factorization 
theorem are derived. Defining the LP and NLP functions, 



Aofa) = m% I dx I d 2 k T ^ kT) 



Ai(ry) = m 2 B J dx J d 2 kx 
respectively, we obtain the form factors, 



cj) + (x,k T ) 



T)m B 



x<t>- (x, kr) 
rj{rjxm B + k\ 



(24) 



F v , A (q 2 ) = — lMv)±^i(v)] 

TUB 



(25) 



Because of kr ~ O(A) in the B meson, Ai(7y) is of 0(A/rrib) relative to Ao(r/) in the large r\ region. Again, only a 
single B meson wave function is relevant in the LP analysis of the B — > decay, consistent with the observation in 
[49|. Compared to Eq. H2l)fl . both (j)± appear in the kr factorization theorem at NLP. 
The decay spectrum is then given, according to Eq. ©, by 



dT 
drj 



48tt 2 



(26) 



Similarly, the NLP terms have cancelled, and only the 0(A 2 /m 2 ) term A 2 (r)) is left. We adopt the models for the B 
meson wave functions in |l3j . whose kj- dependence is coupled to the x dependence through a ^-function, 



1 



4>±(x, Ht) — 4>±(x) — S (fe|. — x(2\ ~ x^m 2 ^) 



(27) 



Using the same input parameters, we obtain the photon energy spectra from fcx factorization theorem in Fig. 0] 
for A = 0.1 and for A = 0.15. These spectra are symmetric in r], and modified only slightly by the higher-power 
correction. Hence, the higher-power correction is under control in fcr factorization theorem compared to that in 
collinear factorization theorem: the power behavior 1 /rj of the spectrum in the small rj region has been smeared into 
r/ln 2 rj. It implies that the perturbative QCD (PQCD) approach based on kr factorization theorem [42l l43l lil. |45| 
has a better convergence at subleading level. 
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FIG. 4: Spectra Spectra in units of GeV 1 from kr factorization with the solid (dashed-dotted)line corresponding to the LP 
contribution for A = 0.1 (A = 0.15), and the dashed (dotted)line to the inclusion of the NLP contribution for A = 0.1 (A = 0.15). 



III. LONG- AND SHORT-DISTANCE CORRECTIONS 

In this section we discuss the long-distance and short-distance corrections to the B — > jlu decay spectrum. For this 
purpose, the form factors are written, in k,T factorization theorem, as 



h 



FvAf) = — AofaHA^fa) +F^ A {q 2 ), (28) 

tub L J 

where Aq 1 "* and Fy 1 ^ denote the 0{a s ) and long-distance correction to the leading result, respectively. We shall 
estimate the latter by considering the B — > V — > 7 transition. This correction is certainly significant in the small 77 
(large q 2 ) region, where the internal quark becomes soft, and easily form a resonance with the spectator quark. Hence, 
it could break the QCD factorization of the form factors Fy,A a t small n. At large 77, the long-distance contribution 
may be suppressed by the values of the B — > V transition form factors 1151 ]. 
The long-distance amplitude is written as 0] 

\{ 1 {P 2 ,e T )\u^{l- 1 ,)b\B{P l )) = (7 (P2,e T ) \Jem\V (P 2 ,e T )) p , _ ^r—^ 

x(V(P 2 ,e T )\u lfl (l- 75 )b\B(P 1 )) , (29) 

with the vector mesons V = p, uj, ■ ■ ■, and their masses my and widths Ty. Take the B meson transition into a 
transversely polarized p meson as an example, for which the first matrix element on the right-hand side of Eq. 129(1 
gives 

< 7 (F 2 , e T )\ J: m \p(P 2 , e T )) = - % ~m p f p ^ , (30) 
f p being the p meson decay constant. The second matrix element is decomposed into 

(p(P 2 ,e T )\uip(l-l 5 )b\B(P 1 )) = - 2V{ f ] e^patfPfPZ -i(m B +m p )A 1 (q 2 )e* Ttl , (31) 

TUB + TTlp 

with the B — > p form factors V(q 2 ) and Ai{q 2 ). Combining Eqs. l(3U|) and l(3T]l. we extract from Eq. 

F^(q 2 ) = /P T ^ V(q 2 ) , 
m p — il p m B + m p 

FT{q 2 ) = fp . r (mB+Wp) A 1 (g 2 ) . (32) 
m p - i\ p r\m B 
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For the long-distance contribution through the B — > uj transition, we have the similar expressions to Eq. Il.'i2[| . but 
with the charge factor 1/2 in Eq. l|30[l being replaced by 1/6, and the appropriate replacement of the vector meson 
mass and of the decay constant. The B — > ip transitions do not contribute in this case. 
For the p and lo mesons, we employ the inputs [4r| 

m p = 0.771 GeV , T p /m p = 0.21 , f p = 0.217 GeV , 

m u = 0.783 GeV , T w /m w w , f u = 0.195 GeV . (33) 

For the B — ► p, cj form factors, we adopt the models derived from the light-front QCD 47j, which have been 
parameterized as 



F{l 2 ) = z ? 2/ ff?,, 2/ 2 , 2 , (34) 

1 — a{q z /m B ) + b{q z /m B y 



with the constants, 



V(q 2 ) : F(0) = 0.27, a = 1.84 , 6=1.28, 

Ai(g 2 ) : F(0) = 0.22, a = 0.95 , 6 = 0.21. (35) 



We restrict the above formalism in the region, 



v > i _ <k™L = 0.275 , (36) 
mi 



Imax being the maximal value of q 2 in the B — > w transition, in which Eq. H34|l holds. The long-distance contribution 
increases Fy t A by about 30 ~ 50% for A = 0.1 ~ 0.15 at large 77, consistent with the observation in [ili I25I Esj. Its 
effect to the decay spectrum is quite important, especially for 77 < 0.8, as shown in Fig. [5j 

The B — > p,u> transition form factors at large recoil could be regarded as an 0(ct s ) object This observation 
hints that we should attempt to take into account the NLO short-distance correction to Fy.A- The NLO correction 
to the B — > jlv decay has been computed by several groups [23, [H, [24| m collinear factorization theorem (SCET or 
QCDF). However, we need the result from kx factorization theorem (with the parton transverse momenta fcy being 
included), which is quoted from [2(ij |: 



(!) g s (2E 7 ) 2 f f 2 <j) + (x,k T ) 

A ° {v) = — —° Fm B J dx J d fc ^ ^|Tfc| 



, 2 77 5 , 77 47T 2 o / k% \ , / fc2, 

In 2 - - - In - + — - In 2 1 + — + 2ttt: In 1 1 T 



X 



2 x 3 V 2fc+ 2 / V 2fc+ 2 



(37) 



The weaker evolution of fs will be neglected for simplicity. Due to the large negative double logarithm, the NLO 
correction to the form factors Fv,a is destructive, and about 30% of the leading result for both A = 0.1 and A = 0.15 
at large 77. The resummation of this double logarithm to all orders has been discussed in [2(1 1221 1231 124 128| . 

We emphasize that the NLO hard kernel depends on a factorization scheme, in which the B meson wave function 
is defined [23^| . Therefore, it is not very legitimate to adopt an expression straightforwardly from some other works in 
the literature. The calculation of the NLO hard kernel for the B — > jlu decay in the factorization scheme specified 
in [3{| is in progress, which will be published elsewhere. The NLO correction in SCET has been further factorized 
into a function characterized by the scale nib, and another by \/ m&A. As stated in [23l |. this further factorization 
is not numerically essential for nib ~ 5 GeV. On the other hand, the model-dependent estimate of the long-distance 
contribution also suffers large uncertainty. Hence, we just intend to point out the potential strong cancellation between 
the long-distance and short-distance corrections in this mode. As shown in Fig. |3J after combining both subleading 
contributions, the net effect has been greatly reduced. Especially, for the shape parameter A = 0.1, the cancellation is 
almost exact for 77 > 0.8. We conclude that the leading result in the large 77 region is stable under these corrections. 

Using the lifetime of a charged B meson t b ± — 1.674 x 10~ 12 s and considering only the leading contribution, we 
obtain the branching ratios for A = 0.15 ~ 0.1, 

B(B -> -/lis) = (1.8 - 4.8) x 10~ 6 , (38) 

from Eq. I|26() in fc-r factorization theorem (PQCD), with only the 0(A 2 /m 2 ) and 0(a 2 ) uncertainty. The values 
in Eq. I|38|l are more or less consistent with other estimates in the literature: a model-dependent evaluation of the 
structure-dependent photon emission contribution gave the branching ratio 10~ 7 ~ 10 -6 01- Using the B meson 
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FIG. 5: Spectra in units of GeV -1 for A = 0.1 and A = 0.15 with the solid lines corresponding to the LP contribution only, the 
dashed-dotted lines to the inclusion of long-distance contribution, the dashed lines to the inclusion of the NLO correction, and 
the dotted lines to the inclusion of both the long-distance and NLO contributions. 



bound-state wave function from a Salpeter equation, 0.9 x 10~ 6 has been obtained [Tt], |2l|. Both a simple non- 
relativistic model and light-front QCD led to 3.5 x 10~ 6 0,0]. Light-cone sum rules and the pole-model calculation 
gave 2 x 10 -6 and 2.26 x 10 -6 (2||, respectively. At last, the experimental upper bound at 90% confidence level 
was oO 

B(B -► jlv) < 2.0 x 10~ 6 . (39) 



IV. CONCLUSION 



In this paper we have studied the B — ► decay in the PQCD approach based on kx factorization theorem. This 
formalism is well-defined at subleading level, since the two-parton LP B meson wave functions remain normalizable 
even after including the evolution effect. Note that the QCDF approach based on collinear factorization theorem 
fails at NLP. We have shown that the 0(l/mb) corrections from the heavy-quark expansion can be absorbed into the 
LP B meson wave functions redefine by the nonlocal matrix element in Eq. i|14|) . The NLP contributions from the 
hard kernels to the decay spectrum cancel. The three-parton B meson wave functions turn out to be suppressed by 
X/m\ in this special mode. The constructive long-distance contribution almost cancels the destructive NLO radiative 
correction for both the form factors Fy and Fa- The B meson wave function 0+ can then be extracted from the 
observed B — ► jlv decay spectrum using the leading formalism, which suffers only the 0(l/m|) and 0{a s ) uncertainly. 
We conclude that the B — * jlu decay is the cleanest mode for determining this important nonperturbative input for 
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the perturbation theories of exclusive B meson decays. The determination can be refined by including the evolution 
and resummation effects into the factorization formulas poll22ll23.l23, l28| . 
Measuring the B — ► spectrum in the lepton and photon energies [20j, 

j?T_ = °&F\V*F< (1 _ ^ {[F , (g2) + F 2 {q 2 )m _ y){1 _ y _ v) + ft 

-2F v (q 2 )F A (q 2 ) V (2~2y~r,)} , (40) 

with the lepton energy fraction y — 2Ei/niB, 1 — ^7 < J/ < 1, we can extract the information of the form factors Fy 
and Fa separately. It is then possible to fix the two two-parton B meson wave functions <p± simultaneously from 
Eq. (|25[) . At this NLP level, the three-parton wave functions are still absent following the reasoning in Sec. II. B. 
The long-distance contribution and the NLO corrections also cancel each other as indicated in Eq. (|28[) . With the 
B — > jlu branching ratio around 10~ 6 , the above experimental determination is possible. 

We thank I. Bigi, Y. Kwon, Z. Ligeti, M. Neubert, T. Onogi, and A.I. Sanda for helpful discussions. HNL acknowl- 
edges the hospitality of Nagoya University during his visit, where this work was initiated. This work was supported 
by the National Science Council of R.O.C. under Grant No. NSC-93-2112-M-001-014. 



APPENDIX A: THREE-PARTON CONTRIBUTION 



We start with Eq. (1.3) in Ref. |5l| : 

G^)= [d^-^^i 9 AM^L, (Ai) 

J 2-k z [z — wY 2n^w^ 
which describes the interaction of a quark with a gluon. In momentum space the above expression becomes 

& [Z) -J (2tt) 4 J (27r)4 e (fc 2 +0 27 I 2 9 ° { 2 > ' [ ' 

where I {k^) is the momentum carried by the incoming quark (gluon). The Feynman parametrization gives 

g<,, <*> - -h /^' (r+ *y . (A3> 

where the variable change / + uk2 — > I, u = 1 — w, has been applied. 

In the case we are considering, the gluon momentum k% is of 0(A), since the B meson is dominated by soft dynamics. 
We expand the above expression up to 0(&2): 

= - /«« ^l^ + ^-^|<A«, 

u fi+~f a j/~ u ]h/ a fi + 

T 3 





f d 4 k 2 




' (2tt) 4 


dH M 


' vr v. 


(27T)4 e \ 


(i 2 ) 2 11 



du / 7TZu e { 77^T*S 'A a {uz) + — ^ T igd A a (uz)n p _ \ . (A4) 



The first term on the right-hand side of Eq. l|A4fl . contributing to a phase factor pSl). will be dropped. For convenience, 
we work in the light-cone gauge A + = 0, in which the second and third terms are rewritten as 

G^iz) = i J dugG aP {uz)ni J ^' "^^^^ . (A5) 

It is clear that the field strength gG a p(uz)n^_ can be factored together with the rescaled b quark field h and the light 
quark field u into the nonlocal matrix element in Eq. I|15|l. The integrand depending on I is then identified as the 
hard kernel in momentum space for the three-parton contribution. Employing Eq. 1A5JI for Fig. [21 and substituting 
P2 — k\ — uki for I, we obtain Eq. IjRjfl . 
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